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ABSTRACT 

We investigate the evolution of 5 granular-scale magnetic flux cancellations 
just outside the moat region of a sunspot by using accurate spectropolarimet- 
ric measurements and G-band images with the Solar Optical Telescope aboard 
Hinode. The opposite polarity magnetic elements approach a junction of the in- 
tergranular lanes and then they collide with each other there. The intergranular 
junction has strong red shifts, darker intensities than the regular intergranular 
lanes, and surface converging flows. This clearly confirms that the converging 
and downward convective motions are essential for the approaching process of 
the opposite-polarity magnetic elements. However, motion of the approaching 
magnetic elements does not always match with their surrounding surface flow 
patterns in our observations. This suggests that, in addition to the surface flows, 
subsurface downward convective motions and subsurface magnetic connectivities 
are important for understanding the approach and collision of the opposite polar- 
ity elements observed in the photosphere. We find that the horizontal magnetic 
field appears between the canceling opposite polarity elements in only one event. 
The horizontal fields are observed along the intergranular lanes with Doppler red 
shifts. This cancellation is most probably a result of the submergence (retraction) 
of low-lying photospheric magnetic flux. In the other 4 events, the horizontal field 
is not observed between the opposite polarity elements at any time when they 
approach and cancel each other. These approaching magnetic elements are more 
concentrated rather than gradually diffused, and they have nearly vertical fields 
even while they are in contact each other. We thus infer that the actual flux 
cancellation is highly time dependent events at scales less than a pixel of Hinode 
SOT (about 200 km) near the solar surface. 

Subject headings: Sun: granulation — Sun: magnetic fields — Sun: photosphere 
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of sizes, ranging from the kilo-gauss fibrils with diameters of a few hundred kilo-meters, 
through the mesoscale pores, to sunspots of various sizes. In the turbulent convection zone 
immediately below the quiet-Sun photosphere, the magnetic field is dominated by the fluid. 
Therefore, it is expected that the field distribution is organized into the convective cells. 
This field distribution is necessarily complex with flux connectivities of various scales, across 
individual and multiple cells. Magnetic concentrations that have vertical and strong (kilo- 
gauss) field are observed along the boundaries of the convective cells. A most prominent 
feature in the quiet Sun is the network magnetic field that partially outlines supergranular 
cells. It is believed that the network magnetic field is formed by the advection of internetwork 
fields via the supergranular flows. Magnetic elements in internetwork areas tend to move 
toward the nearest network concentration at a speed of about 0.2 km s~^ (Ide Wiin et al. 



2003 



2008 ) , and their rms velo city is about 1.5 km s ^ in internetwork areas (INisenson et al. 



de Wijn et al.ll2008l ). Velocities of 0.2 km s ^ and 1.5 km s ^ are similar to a typical 



speed of supergranular flows and granular flows, respectively. 

Down at the scales of a few hundred kilo-meters, equal to a few photospheric density 
scale heights, magnetic elements of both polarities form to often merge and collide with each 
other in the photospheric fluid. The net magnetic flux increases or decreases, depending on 
whether like-polarity or opposite-polarity magnetic elements have merged. In this paper, 
we report an interesting set of "magnetic flux cancellation" events of this kind observed 
at the high spa tial resolution of t he Solar Optical Telescope (SOT; iTsuneta et al.l |2008[ ) 
aboard Hinode f lKosugi et al.l 120071 ). The magnetic flux cancellation is a descriptive term to 
indicate a mutual flux loss due to the apparent collision of the opposite-polarity magnetic 
elements (IMartin. Livi. fc Wangl 119851 ). This flux cancellation is essential to the process of 
replacement of old magnetic flux w i th newly emergi ng flux in the quiet Sun on a timescale 
of a few days (jSchrijver et al.l Il997l : iHagenaarl 120011). and also to the process of removal of 
sunspot magnetic flux from the photosphere (iKubo et al.ll2008l ). 



Vario us possible processes have been p roposed to explain the observed flux cancella- 
Ryutova et al.l |2003[ ) . involving submergence (retract) of fi-shaped 
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tion (e.g. 

loops or emergence of U-shaped loops across the photosphere. In both cases, the canceling 
opposite-polarity magnetic elements correspond to the two intersections of such loops with 
the photospheric layer. The opposite-polarity magnetic elements disappear w hen the top o f 
a submerging f2-loop has passed through the photospheric layer (see Fig. 2 in IZwaanlll987l ). 
Alternatively, these elements disappear when the bottom of a rising U-loop has passed clear 
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through the photospheric layer. When magnetic field lines have emerged into the chro- 
mosphere and corona, they can hardly submerge back below the photosphere because of 
magnetic buoyancy. Magnetic reconnection taking place within several scale heights above 
the solar surface is probably needed to create lo w-lying O-loop s whose magnetic tension force 
can then overcome the magnetic buoyancy force (jParkeiill975l ). In the photospheric magnetic 
reconnection cases, reconnection should be mo st efficient around the temperature minimum 
region: about 600 km above the solar surface (ILitvinenkd Il999l : iTakeuchi fc Shibatall200ll ). 
In contrast, both magnetic tension and buoyancy forces are directed upward in the case of 
a U-loop rising through the photosphere. Magnetic reconnection is therefore not crucial for 
the emerging U-loops. 

An important observable signature for understanding fiux cancellation is the motion 
of the horizontal magnetic field connecting the canceling magnetic elements. Horizontal 
magnetic fields have been observed between the opposite-polarity magnetic elemen ts during 
the cancellations of moving magnetic features around a sunspot (jChae et al.l 120041 ). Similar 
horizontal fields have a lso been observed in events of cancellations of pores and sunspots 
( iKubo fc Shimizull2007l ). A fiux cancellation witho ut increase of the horizonta l field has also 
been reported for the moving magnetic features (IBellot Rubio fc Beckl l2005l ). Knowledge 
of the full vector field permits one to determine whether the field geometry has fi-loop 
topology or U-loop topology, but in the case of cancellation of small, isolated fiux elements, 
such a determination is often compromised by usual difficult y of resolving t he 18 0° azimuth 
ambiguity. Regarding the motions at such a cancellation site, iHarvey et al.l ( 119991 ) show that 
the magnetic fiux disappears in the chromosphere before it does in the photosphere for at 
least about half of the cancellation events. They suggest that magnetic flux is submerging 
i n most, if not al l, of the cancellation site s. On the other harid, bo th Doppler red shift 
(IChae et al.l 120041 ) and Doppler blue shift (lYurchyshyn fc Wangll200ll ) are reported in the 
cancellation sites. The center-to-limb variations of the Doppler velocities at the polarity 
inversion lines in the cancellation sites suggest that the observed velocities with the spatial 
resolution of about 1" mainly show the gas flows al ong the horizontal fields rather than actual 
emerging or submerging motions of the field lines ( IKubo fc Shimizd 120071 ). Observationally, 
the physical process of the magnetic fiux cancellation is still not well understood. 

Recent observations with the high spatial resolution by Hinode SOT show that down- 
ward motions (^r ed-shifts) are continuously observed during the fiux cancellation process 
( llida et al.l 120101 ). From the Stokes-V signals far from the disk center, they also suggest 
that the canceling opposite polarity elements tend to have an f2-shaped configuration rather 
than a U-shaped configuration. In this paper, we investigate the detailed evolution of can- 
celing magnetic elements and their surrounding convective motions near the disk center by 
using vector magnetic fields and velocity fields observed with the SOT. We are particularly 
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interested in the evolution that precedes the cancellation in order to investigate why the 
opposite-polarity magnetic elements approach and collide with each other. For a full under- 
standing of the phenomenon, this dynamical development is as important as the flux removal 
process at the cancellation site. 



OBSERVATIONS AND ANALYSIS 



Our target flux-cancellation events occurred just outside an active region NOAA 10944, 
as shown by the solid boxes in Figure [H We selected a data set simultaneously taken by 
the SOT spectropolarimeter (SP) and filtergraph (FG) in their full spatial resolution modes. 
The SP repeatedly scanned the same area with a field of view of 10" x 82" (the dotted box 
in Fig. [I]) for 4.5 hr from 06:35:32 on 2007 March 2. The SP measured Stokes I, Q, U, and V 
profiles across Fe I 630.1 nm and 630.2 nm lines. A width of a slit was 0".15 and a pixel sampling 
along the slit was 0".16. One scan took 5.5 minutes with an integration time of 4.8 s at each 
slit position. We estimated total circular polarization (Ctot) and total linear polarization {Ltot) as 
follows: 
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where Aq is the center of Fe I 630.2 nm line in each pixel, and 1^ is the local continuum intensity. 
The continuum intensity is average d over the Stokes I p rofile from Aq + 43.2 pm to Aq + 64.8 pm. 
In the weak field approximations (jJefferies et al.lll989l ). Ctot is proportional to the longitudinal 
magnetic flux density {fBi) and Ltot has the relation to the transverse flux density as Ltot ~ 
where Bl, Bt, f are the longitudinal field strength, transverse field strength, and filling factor, 
respectively. The observed area is located not far from the disk center (S06W21). We can consider 
the longitudinal and transverse fields to be the vertical and horizontal fields with respect to the 
solar surface, respectively. In addition to the magnetic field, we compute the Doppler velocity with 
respect to the average of Doppler velocities in the quiet area. 

During the repeated scans by the SP, the FG took a time series of G-band images with a pixel 
scale of O". 05448. The cadence is 1 minute, and the field of view is 56" x 112". After applying a sub- 
sonic filter for the time series of G-ba nd images, we estima.ted gr anular flow patterns with the local 
correlation tracking technique (LOT: [November &: Simonlll988l ). The apodization window of the 
LCT had a Gaussian shape with the FWHM of 32 pixels (l".7). The LOT flow maps are calculated 
from the series of G-band images. The cadence of the LCT flow maps is therefore 1 minute. Drift 
of the image, arising from the correlation tracker ([Shimizu et al.ll2008l ). was removed by aligning 
the sunspot centers in the G-band images at different times. After removing the image drift, both 
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G-band images and flow maps were averaged over the period corresponding to the duration of an SP 
map (5.5 minutes). The SP maps were ahgned to the averaged G-band images at the time closest 
to the midpoint of the SP maps. The ahgnment was performed by the image cross-correlation 
between the averaged G-band images and the SP continuum intensity maps. Thus, the image drift 
in the time series of the SP maps was removed as well as the G-band images. 



3. RESULTS 

Here we present the evolution of flux cancellation events in Region A and Region B of Figure [1] 
We shall henceforth use simplifled terminology of calling the positive-polarity and negative-polarity 
magnetic elements simply as the positive and negative elements, respectively. 



3.1. Region A 

The flrst event is an elemental flux cancellation in Region A. Figure [2] shows the developments 
of this region prior to the start of the cancellation event. Fuzzy positive (white) and negative (black) 
elements appear within the dashed circle in the map of the total circular polarization {Ctot}- In 
this circle, an increase of the total linear polarization (Ltot) is also observed. These results indicate 
that the opposite-polarity elements emerge into the photosphere as a pair whereas the canceling 
magnetic elements in the dashed box do not originally emerge as a pair. The emerged negative 
elements merge together to consolidate into a stronger negative element that subsequently approach 
and mutually cancel with a positive element. Even when the negative element is close to the positive 
element, there is only a slight increase of Ltot in a part of the region between them. 

Figure[3]a shows that, like the negative canceling element, smaller positive elements around the 
cancellation site first converge to form into a prominent positive element that then moves toward 
the negative element. We trace the center of the positive element. The center (the cross symbol 
in Fig. 3a) is defined as the average of the position (r) weighted by Ctot- Yl^tot^i/ Yl^tof This 
center moves at about 1.4 km s~^ during 11 minutes. Note tha t the motion of 1.4 km s~^ is similar 



to th e rms velocity of magnetic elements in internetwork areas (jNisenson et al.ll2003l : Ide Wijn et al 



20081 ). The positive element collides with the negative magnetic element at a time between the 
third and fourth frames. The negative element disappears from the photosphere by 10 minutes 
after the start of the cancellation, whereas the positive element does not completely disappear in 
this event. Even while the net flux of both negative and positive elements is decreasing, Ctot of 
each pixel in the colliding magnetic elements does not decrease. Furthermore, both Ctot and Ltot 
slightly increase at the center of the positive element at the fourth frame in Figure [3l This increase 
arises from either the increase of the fllling factor or the held strength within the positive element. 
These results suggest that the colliding magnetic elements have nearly vertical magnetic flelds 
until these magnetic elements disappear form the photosphere. On the other hand, no increase of 
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Ltot is observed between the canceling opposite polarity elements during the whole flux cancellation 
process. This means that the horizontal magnetic fields, which are expected in both the submerging 
il-loop model and the emerging U-loop model, do not appear at the cancellation site. Figure H] 
confirms that neither Stokes Q nor U signal increases during the fiux cancellation. 

The cancellation site is located at a junction of the intergranular lanes (the upward arrow in 
Fig. [3lfi). The intergranular junction is darker in G-band intensity than its nearby intergranular 
lanes. Intergranular lanes have generally red shifts, and a larger red shift with about 1.5 km 
s~^ is observed at the intergranular junction just before the cancellation (Fig. [3]c). Interestingly, 
such a larger red shift is not observed there during the cancellation. The left panels of Figure H] 
also show that the red shift of Stokes I profile is only seen in the top panel taken just before 
the cancellation. This means that the observed large red shift at the cancellation site does not 
result from this fiux cancellation process. Instead, we may infer that the opposite polarity elements 
approach the intergranular junction with the large red shift. The intergranular junction tends to 
have a converging flow rather than a diverging flow, but the flow pattern is not clear in this event. 
The size of the intergranular junction is too small to examine the flows derived by the LOT method 
because the LOT window is 1".7. One interesting result is that the negative element moves toward 
such a small, less clear converging area although a large, strong converging area with the darker 
intensity is also located near it, as shown by the dashed circle in Figure [3]e. 



3.2. Region B 

Figure [5] shows a negative element in Region B that evolves into three branches. Each of these 
three branches elongates toward an area with a large red shift (the arrows in Fig. [5|). Magnetic 
flux cancellations subsequently occur at the tip of each branch (the dashed circles in Fig. [5]). An 
increase of Ltot signals indicating horizontal magnetic fields is observed only in the cancellation event 
enclosed by the final dashed circle. The other 3 cancellation events do not accompany the horizontal 
magnetic fields, and their properties are similar to the fiux cancellation observed in Region A. In 
the event with the increase of Ltot, the horizontal fields are located along the intergranular lane 
with the Doppler red shifts. The sizes of the magnetic elements begin to decrease slightly in the 
Ctot maps just before their collision. This event is different from the other cancellation events: the 
horizontal field {Ltot signals) is already present as a result of the fiux emergence well before the 
event of cancellation sets in, as identified by the dashed box at the beginning of the observation 
sequence on the top of Figure [5]a. This fiux emergence has a possibility of distinguishing this event 
from the other four cancellation events. However, it is unclear whether these horizontal fields are 
the same as the horizontal fields observed during the cancellation or not, because the Ltot signals 
in the dashed box disappear briefiy before the opposite polarity elements approach each other. 

The circles in Figure[6]show that the areas with strong converging fiows mostly have a darker G- 
band intensity in the intergranular lanes. The size of these converging fiow areas is larger than that 



-7- 



of the cancellation event in Region A, but still smaller than a travel distance of the approaching 
magnetic elements. The areas with strong converging flows correspond to the junctions of the 
intergranular lanes, and arc observed to have downward Doppler motions. The point to emphasize 
is that strong converging areas do not necessarily have cancellation events, whereas the cancellation 
sites tend to be located beside or at the strong converging areas, as suggested by our observations. 

4. DISCUSSION 

We have observed the detailed evolution of magnetic fields and velocity fields for 5 cancellation 
events at the granular scales just outside the sunspot moat region. Prom the observations at high 
spatial resolutions, we have clearly confirmed that opposite-polarity magnetic elements mutually 
cancel by moving toward the junctions of intergranular lanes. The intergranular junctions arc 
characterized with darker intensities than their nearby intergranular lanes, strong red-shifts, and 
surface converging flows. Our new finding is that no horizontal field (no increase of Ltot) appears 
between the canceling opposite polarity in 4 of the 5 flux cancellation sites. In the events without 
the horizontal field, the canceling opposite-polarity magnetic elements have nearly vertical fields 
even while their net magnetic flux decreases. Here, we discuss the physical issues posed by these 
cancellation events, for flux removal and flux-tube collisions as dynamical processes. 

4.1. Cancellation without Appearance of Horizontal Fields 
4.1.1. Flux Removal from the Photosphere 

The submergence of an $7-loop is dynamically quite different from the emergent rise of a U- 
loop. Nevertheless, the photospheric convergence of loop footpoints in both processes produces a 
cancellation of opposite elements with the same magnetic signature. In each case, the cancellation 
of opposite vertical fields is accompanied with an increase followed by a decrease in the horizontal 
field. This magnetic signature is not found in our observations. The flux cancellation process 
in our observations is not spatially resolved even with the Hinode SOT. In other words, flux is 
removed from the photosphere at sizes smaller than the ~200 km spatial resolution of the SOT. The 
observed Doppler velocities at the cancellation sites are probably related to neither submergence nor 
emergence of magnetic field lines. These velocities are just the expected downdrafts of the convective 
motions at the intergranular lanes which are present even when flux cancellation is not occurring. 
We believe that the absence of the horizontal field in the cancellation sites is not due to cither a 
low cadence or an insufficient signal-to-noisc ratio (S/N). The decrease of the canceling (colliding) 
magnetic elements is usually observed in two or more frames. There is a low possibility that the 
horizontal field appears only during the period when the slit is located outside the cancellation 
sites. Moreover, small-scale horizontal magnetic fields are observed outside our cancellation sites 
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in th i s data set. Such horizontal fields typically ha ve an order of hecto-gauss (jOrozco Suarez et al 



20071 : iLites et al.l 120081 : llshikawa &: Tsunetal 12003 ) . Horizontal magnetic fields having sizes larger 



than the spatial resolution must be detected at the cancellation site unless such horizontal fields 
suddenly become one order of magnitude weaker than the vertical fields in the colliding magnetic 
elements. 



4.1.2. Approach and Collision of Opposite-polarity Magnetic Elements 

Another important issue is why a magnetic polarity element would move in such a manner as to 
meet with another element of the opposite polarity on the vast solar surface. One possibility is the 
chance encounter of these elements advected by either granular fiows or supergranular flows. The 
magnetic elements advected by these flows are likely to stay in the intergranular lanes, especially 
the boundary of supergranular cells. Flux cancellation events also prefer to occur there. This 
is consistent with our observations that magnetic fiux cancellations occur in the intergranular 
junctions having the strong converging flows. However, one remaining issue is that motion of the 
approaching opposite polarity elements is not necessarily consistent with the surrounding surface 
flow patterns. The approaching opposite-polarity elements travel a longer distance than the size of 
the strong converging areas, and do not always move toward the nearest, stronger converging area. 
Moreover, it is difficult to explain such motion of magnetic elements from only the advection by 
supergranular flows, because these magnetic elements move at the speed similar to usual granular 
flows, which is faster than a typical speed of supergranular flows. One possibility is that we have 
missed high-speed, systematic flows along the intergranular lanes because of the insufficient spatial 
resolution in our LOT velocity maps. Although such flows are not yet reported, supers onic granular 



horiz ontal flows recently detected by Doppler measurements far from the disk center ([Bellot Rubio 



20091 ) may be able to produce the high-speed fiows along the intergranular lanes. Another possibility 
is that other forces (e.g., an intrinsic Lorentz force) in addition to the force driven by the surface 
flows are needed to explain the motion of the approaching magnetic elements in the photosphere. 

The subsurface advection of magnetic fields is as important as the surface advection for think- 
ing about the observed motion of photospheric magnetic elements. In particular, if a subsurface 
field connects the opposite-polarity magnetic elements as a U-loop, our observations suggest the 
following interesting possibility. A cooler material sinking to the bottom of the U-loop would not 
only prevent the loop from rising, but may even force the loop to sink with the converging and 
downward convective fiows, as sketched in Figure [71a. The cooler material can drag the subsurface 
field lines into a deeper layer via a downward convective flow because the plasma-/3 is higher than 
unity below the photosphere, except for the strong magnetic fields associated with sunspots. This is 
similar to the idea of downward fiux pumping by the turbulent granular convection around sunspots 



(jThomas et al.ll2002l ). As a result of the forced submergence of the subsurface U-loop connecting 
the opposite polarity elements, these magnetic elements are driven by their Lorentz force to move 
toward the area with the large downward motion and then collide with each other (the middle 
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panel of Fig. [71a). More detailed behaviors of flux tubes with the submergence of the subsurface 
U-loop can be found in Appendix A. Such a process is effective not only for the collision of the 
opposite polarities but also for the merging of the same polarities that are observed before the 
cancellation. Magnetic reconnection just above or below the solar surface would be needed for the 
disappearance of the colliding opposite polarity magnetic elements from the photosphere (the right 
panel of Fig. [7]a), because the bottom of the subsurface field lines dragged by the downdrafts of 
cooler materials does not easily emerge into the photosphere. The reconnection site is still unknown 
from our observations, but the unresolved fine-scale flux removal process at such a site suggests 
magnetic reconnection, if any, is close to the photospheric surface. We of course have no obser- 
vational evidence that a subsurface U-loop connects the opposite-polarity elements. Nevertheless, 
the events observed in this study can be explained by such a subsurface process. 



4.2. Cancellation with Appearance of Horizontal Fields 



The horizontal magnetic fields between the canceling opposite polarity elements have been 
observed along the intergranular lanes characterized with Dopple r red shifts. The horizontal fields 



with t h e red shifts at the cancellation sites are already reported bv lChae et al.l (j2004l ): ICheung et al 



(|2008l ) ; llida et al.l (|20ld ) . This cancellation is most probably a result of the submergence (retrac- 
tion) of low- lying photospheric field lines along the intergranular lanes (Fig. [7]&). Such low- lying 
photospheric field lines probably are formed by magnetic reconnection in the photosphere (or the 
bottom of the chromosphere) or by a failure of the emergence into the upper atmosphere. Recent 
three-dimensional (3D) radiative MHD simulations show th e granular-scale flux cancellation due 
to a retraction of the 0-loop within an emerging flux region (jCheung et al.l 120081 ). In t he previous 



studi es, the flux cancellations with the h orizontal fields are r eported in the moat region (jChae et al 



2004 ). within the emerging flux regio n (jCheung et al.ll2008l ). or around the center of complicated 
active regions (jKubo &: Shimizull2007l ). These regions basically contain many horizontal fields even 
if the horizontal field is not observed just before the cancellation. In our cancellation event with 
the horizontal field, small fiux emergence accompanied by the appearance of horizontal fields is 
also observed well before the start of the event. The observing products during the cancellation 
may depend on the magnetic field configuration already formed before the approaching and can- 
celing process. On the other hand, the opposite polarity elements approach the region with the 
converging fiows, the darker intensity, and the large red shifts as in the case of cancellations without 
the observation of the horizontal fields. Therefore, the converging and downward motions are also 
important to bring one magnetic polarity element to another polarity element in this event. 
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5. CONCLUSION 

We have presented an observational study of flux cancellation events on the photosphere that 
are sufficiently resolved by the Hinode SOT to show that they are characterized by the absence 
of a horizontal field during the cancellation process. These events are interesting because in the 
usual idea of the submergence of a low-lying O-loop or the buoyant rise of a U-loop, the appearance 
of a horizontal field is the observational signature of the loop top (or bottom) passing across 
the photosphere. Such flux cancellations appear to be more common than the cancellation with 
the appearance of the horizontal fields. Although the nature of the flux removal process at the 
cancellation site is an open question for cancellations without the appearance of the horizontal 
fields, our study shows that it takes place in local areas at scales less than the ~200km resolution of 
the SOT and close to the solar surface. The distinction between the cancellations with and without 
the appearance of the horizontal fields might arise from whether the canceling opposite polarity 
elements have emerged into the photosphere as a pair or not. However, we have investigated only 
5 events that have the size less than a few arcseconds just outside the moat region of a sunspot. 
We need to investigate what really determines the appearance of the horizontal fields between 
the canceling magnetic elements by using more events. In particular, the origin of the canceling 
magnetic elements and their surrounding magnetic field configuration (formation of pre-existing 
horizontal fields) may be important. 

We have confirmed the converging and downward convective flows are essential for the ap- 
proaching and canceling process of the opposite-polarity magnetic elements. The flux cancellations 
are observed at the intcrgranular junction characterized by the strong converging and downward 
flows. However, the approaching opposite-polarity elements seem not to always follow the sur- 
rounding surface flow patterns, at least in our observations. We are accustomed to thinking of 
the connectivities of the surface magnetic fields that we can observe in the solar atmosphere. Our 
observational study suggests that it is important to also consider the connectivities of the surface 
fields that occur below the photosphere. Lanes of downdrafts of cool fluids entraining magnetic 
fluxes that thread across these lanes below the photosphere must be a common occurrence. Such 
a process is a simple explanation for a pair of opposite-polarity elements to appear to seek each 
other on the photosphere. Information of subsurface convective flows would be extremely helpful 
for better understanding of the magnetic flux cancellation. 

We are grateful to T. Yokoyama, Y. Katsukawa, A. de Wijn, Y. lida, and Y. Fan for useful 
discussions and comments on this paper. Hinode is a Japanese mission developed and launched 
by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) as international 
partners. It is operated by these agencies in cooperation with ESA and NSC (Norway). The FPP 
project at LMSAL and HAO is supported by NASA contract NNM07AA01C. 
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A. Discrete Magnetostatic Flux Tubes 

The dynamical process sketched in Figure[7]is due to an interplay among pressure, gravitational 
and Lorentz forces in the neighborhood of a subphotospheric convective downflow. A proper study 
of this 3D time-dependent process requires numerical MHD simulation. On the other hand, some 
physical insight into that interplay can be seen in 3D static solutions of discrete magnetic flux tubes 
in a stratified atmosphere. Generally, such static solutions also require numerical computation but 
the family of analytical solutions taken from Low (1982) serves our purpose here. 

Consider the magnetostatic equilibrium equations: 

— (V X B) X B - Vp - p^i = 0, (Al) 
An 

V • B = 0, (A2) 

P = P , (A3 

m 

describing the force balance for a magnetic field B in an isothermal atmosphere at temperature Tq 
and stratified by a uniform gravity with acceleration g the — z-direction, using Cartesian coordi- 
nates. We use the ideal gas law ()A3P where k and m denote the Boltzmann constant and the mean 
molecular weight of the plasma. 

A 3D particular solution of these equations can be constructed for a magnetic field of the form 

B = F((/>,x)V0 X Vx, 



= m^)^0.^,-^j, (A4) 

where -F((/>, x) is an arbitrary function of two variables and 

for some constant k^. Direct substitution of this field into the magnetostatic equations shows that 
the equilibrium pressure must take the form 

= (^0 - -^P^i^^ ^)) exp{-koz), (A6) 

introducing a constant Pq and identifying k^^ = kTo/mg as the hydrostatic scale height which is 
of the order of 300 km at the photosphere at a temperature of about GOOOiT. 

This family of solutions is geometrically quite simple although it is three-dimensionally varying. 
Figure[8]a shows the contours of constant (/>(y, z) on a plane of constant x. The lines of force (LOFs) 
are all geometrically the same on each constant-x plane but the field varies with all three Cartesian 
coordinates through its amplitude function -F((/>, x) which is an arbitrary function to be prescribed. 
Take any explicit function a{(j),x). Then, setting o" = o"o, a constant, generates a magnetic flux 
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surface in 3D space. To construct a flux tube of a finite cross section, a suitable functional form 
of a{(p, x) set to a constant uo describes a flux-tube boundary. Then, prescribing F{(j), x) / 
inside the tube and F{(j), x) = in the rest of the atmosphere completes the construction. The 
atmospheric pressure distribution is then given by Equation (jA6p which can also be expressed in 
the form 

P+ — = Poexp{-koz). (A7) 
ovr 

External to the flux tube, B = so that p = Poexp(—koz), the isothermal pressure of the field- 
free part of the atmosphere. Internal to the flux tube, the atmospheric pressure is reduced and 
compensated by the magnetic pressure so that the total pressure is stratifled in the same manner 
as the external atmospheric pressure. This solution is mathematically analogous to the solution 
describing the equilibrium between a field B = B{x, y)z and fiuid pressure in the absence of gravity. 
Equilibrium in this case is satisfied by requiring the total pressure P = p + S^/Svr to be uniform 
in space. We have complete freedom to prescribe the field distribution and use this requirement to 
obtain the associated equilibrium pressure. 

The static balance of forces may be viewed by writing Equation (lAlh in the form 

^ (B • V) B - V (p + ^VSvr) - pgz = 0, (A8) 

describing the balance among the forces of magnetic tension, total pressure, and gravity, represented 
by the three terms, respectively. Since the total pressure is vertically stratified everywhere, it follows 
that this class of eq uilibrium m agnetic fields is characterized by a tension force that is vertically 



oriented everywhere (jLowlll984l ). Where the atmospheric plasma is threaded with a magnet ic field, 



its p ressure and density are reduced, so that this portion of the atmosphere is buoyant (jParkei 



19791 ). The buoyancy force is just the net sum of the second and third terms in Equation (jASp . and 



it is balanced by the remaining te rm which is the magnetic tension force. These solutions include 



linearly stable equilibrium states (|Lowlll982l ). 



To address the interplay among these magnetostatic forces, let us construct three equilibrium 
states, each corresponding to two fi-shaped flux tubes that are joined below the surface 2 = 
idealized to be the photosphere. These three states are displayed in Figure [8]& shown as vertical 
sections of flux tubes boundary cut by a constant-x plane. For each flux tube, its boundaries 
are drawn in different colors to identify them. Each pair of boundaries is the intersection of the 
flux-tube surface a = uq with the constant-x plane. 

Consider the M-shaped flux tube drawn in black, with boundaries (j){y,z) = 0.4,0.45. This 
flux tube is made up of two identically shaped tubes that join together below the origin where 
it is kinked. This flux tube intersects the photosphere z = at four places, giving two pairs of 
opposite-polarity magnetic elements on either sides of the origin. Since these elements are threaded 
by the same flux tube, the two elements on either side of the origin are of opposite polarities. The 
location of these two elements is controlled by the anchoring of the two parts of the flux tube at the 
kink below the origin. The inflnitely long far arms are free to equilibrate with their surrounding 
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fluid. The flux tube is everywhere in equihbrium except at the kink where an upward Lorentz force 
is assumed to be balanced by some force representing the effect, for example, of the ram pressure 
in a downdraft. The kinked part of the flux tube has a vertical thickness of about 0.05A;o. 

Suppose we forcefully move this kinked part of the flux tube from its original depth at about 
O.lfco below z = to about 0.35ko, keeping the same flux distribution inside the tube and the same 
vertical thickness at the kink. The new equilibrium of the flux tube is shown in red. The entire 
flux tube has sunk with that downward displacement of the kinked part. The opposite-polarity 
elements on the two sides of the origin drift apart to seek their opposite elements further away 
from the origin. This behavior can be understood as follows. The flux tube in this case does 
not extend very high into the atmosphere. The tube footpoint separations are of the order of a 
hydrostatic scale height, and the curvature at the tops of the tube is strong enough to produce a 
tension force comparable to the pressure and buoyancy forces at the temperature of about 6000K 
for the photosphere. When the kinked part of the flux tube is submerged further, an immediate 
effect is to increase that magnetic curvature. This enhances the magnetic tension force over the 
buoyancy force, leading to the submerging displacement for the whole flux tube. 

The flux tube in green shows a different behavior as the result of the buoyancy force. In this 
case, the kinked part of the flux tube in black is submerged with two effects. The submergence 
elongates the vertical thickness of the kinked part of the flux tube while the cross section of the 
tube becomes more narrow. This elongation may be produced by a downward displacement that 
is larger at the bottom than at the top of the kinked part of the tube. The three-dimensionality of 
our magnetostatic solution is essential. If the field does not vary in x, we would have a flux layer 
rather than a fiux tube. The entire atmosphere lying above it is trapped and cannot fall through 
without breaking the layer with variation in x. In our case, we have a true fiux tube with a finite 
cross-sectional area. This means that the atmosphere external to the flux tube can yield and flow 
around a rising flux tube, for example. Moreover, the narrowing of the cross section of the tube 
enhances the magnetic pressure in the tube, which n aturally produces a siphon flow along the flux 



tube to drain fluid down the far arms of the tube (jPikel'Neij|l97ll : lThomaslll988l ). The reduced 
weight of the tops of the tube becomes more buoyant and rise to a greater height. This brings 
about equilibrium with the external atmospheric pressure and a final balance between buoyancy 
and magnetic tension forces. The dynamical process builds up stress that pushes the two far arms 
away from the origin while at the same time brings the two opposite-polarity elements on either 
sides of the origin closer together, in a manner similar to the Hinode observations reported in this 
paper. In our Hinode observations the canceling opposite-polarity elements have fields that return 
to the photosphere at distances much larger than 10 hydrostatic scale heights from the canceling 
sites. We expect the buoyancy force to dominate in the atmospheric portion of the field, keeping 
the tops of these fields up in the chromosphere and above. 

This simple theoretical analysis illustrates the interplay among the magnetic and hydrody- 
namic forces in the submergence of a field connection that happens to be located right at a sub- 
photospheric convective downdraft. We of course need 3D time-dependent MHD simulations to 
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examine the dynamics of this process to get a complete physical picture. 
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Fig. 1. — G-band image taken at 08:00:32 on 2007 March 2. The G-band image is normalized 
to the mean intensity of the quiet area outside the sunspot. The dashed box shows a scan 
area by the Hinode SP. The sohd boxes labeled "A" and "B" are identical to the fields of 
view of Figs. |2] and [5l respectively. 
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Fig. 2. — Time series of (a) the total circular polarization (Ctot) and (b) the total linear 
polarization (Ltot) before the flux cancellation in Region A of Figure [U White (black) 
indicates positive (negative) polarity in panel a. The contours in panel b represent the 
±0.01Ctot levels. The units of vertical and horizontal axes are in arcseconds. The dashed 
circle indicates a flux emerging site. The dashed box in the final frame includes our target 
flux cancellation event, and is identical to the field of view of Fig. [31 
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Fig. 3. — (a) Total circular polarization (Ctot), (b) total linear polarization (Ltot), (c) Doppler 
velocity in units of km s~^, (d) normalized G-band intensity, and (e) flow map during the flux 
cancellation in Region A. The upward arrows in panels a-d indicate the strong red- shifted 
area, and the contours are ±0.01 of Ctot- The cross symbols in panel a represent centers 
of the positive-polarity magnetic element to be canceled. The arrows in panel e show the 
horizontal velocities derived with the local correlation tracking method. The upward arrow 
on the right-hand side of the panel is 1.5 km s~^. The background of panel e is identical to 
that of panel d, and the dashed circle surrounds the area with strong converging flows. 
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Fig. 4. — Stokes profiles of tlie Fe I 630.2 nm line at the pixels represented by squares in panel 
a of Fig. El The profiles in panels a-c are taken at 09:39:39, 09:45:02, and 09:50:30 on 2009 March 
2, respectively. These profiles are normalized by the continuum intensity averaged over the quiet 
area. The vertical dashed line represents the averaged position of the line centers over the map. 
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Fig. 5. — Time series of (a) the total circular polarization {Ctot), (b) the total linear 
polarization (Ltot), and (c) the Doppler velocity. The contours represent the ±0.01Ctot 
levels. The dashed circles and the dashed box indicate the flux cancellation sites and the 
flux emergence site, respectively. The arrows in panel c point to the areas with strong red 
shifts. 
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Fig. 6. — (a) G-band images and (b) flow maps with the local correlation tracking method in 
Region B. The field of view and the contours are identical to those of the third row in FigJS] 
The right-handed arrow at the top of panel b is 1.5 km s~^. The dashed circles encircle the 
areas with strong converging fiows. The background of panel b is identical to that of panel 
a. 
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(b) 



reconnection just above the solar surface 



submerging low- laying photospheric field line 

I 



disappearance of magnetic elements 



converging and downward tlows^ 





Fig. 7. — Schematic illustration for a flux cancellation without (panel a) and with (panel 
b) the observation of horizontal fields. 
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Fig. 8. — (a) Contours of constant z), given by Eq.f lA5|) . on a plane of constant x. (b) 
Three states of the M-shaped flux tube that is composed of two r2-shaped flux tubes joined 
at a kink. 



